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Abstract 

We construct a 4-d lattice Dirac operator D using a systematical 
expansion in terms of simple operators on the lattice. The Ginsparg- 
Wilson equation turns into a system of coupled equations for the expan- 
sion coefficients of D. We solve these equations for a finite parametriza- 
tion of D and find an approximate solution of the Ginsparg- Wilson 
equation. We analyze the spectral properties of our D for various en- 
sembles of quenched SU(3) configurations. Improving the gauge field 
action considerably improves the spectral properties of our D. 
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1 Introduction 



Chiral symmetry in fermionic field theories on the lattice may be imple- 
mented by adding a non-vanishing and in the continuum limit irrelevant 
right hand side of the usual anti-commutator of the Dirac operator D with 
75 . In its simplest form the resulting so-called Ginsparg- Wilson equation |^] 
reads (we set the lattice spacing to 1) 

Dj5 + 75^ = Dj^D . (1.1) 

Based on ( |1.1| ) chirally symmetric fermions can be constructed on the lattice 
(for reviews of recent developments see e.g. [0, ^). 

Currently two types of solutions for (|l.lD are known. Neuberger [^] 
gave an explicit construction, the so-called overlap operator based on earlier 
work, the overlap approach to chiral fermions on the lattice |^. The overlap 
operator is given by 

D = 1 - A (a^aY^ with A = 1 + s - Do . (1.2) 



Here Dq is some decent lattice Dirac operator which is free of doublers. 
Typically the Wilson-Dirac operator is used here but also other choices give 
rise to solutions of (1.1); s is a real number which can, at the cost of an 



additional renormalization, be used to optimize the locality properties of D 
P]. The main challenge when using the overlap operator is the computation 



of the square root in (1^), making simulations with the overlap operator 



considerably more expensive than simulations with the Wilson operator. 



The second solution to (1^) is provided by the fixed point Dirac operator 
which goes back to the perfect action approach to lattice field theories 
10[. So far, the fixed point Dirac operator has been computed only in two 



dimensions [11, |T^, |l^ and it has become clear that its construction in 4-d 
is quite a challenging enterprise. 

Recently H a new line of attack for solving the Ginsparg- Wilson equa- 
tion has been proposed. The basic idea is to systematically expand the 
most general lattice Dirac operator D in a series of simple basis operators 
on the lattice. This expanded D is then inserted into the Ginsparg- Wilson 



equation. The product D^^D on the right hand side of (IT) can be eval 



uated using techniques from the analytical hopping expansion. The result 



is an expansion of both sides of (IT); comparing the factors in front of the 
contributions to this expansion one finds that the Ginsparg- Wilson equa- 
tion is mapped to an equivalent system of coupled quadratic equations for 
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the expansion coefficients of D. When using finitely many terms in the ex- 
pansion of D, the system of equations for the expansion coefficients can be 
solved numerically and the corresponding D is an approximate solution of 
the Ginsparg- Wilson equation. By including more terms in the expansion 
of D the violation of can be made arbitrarily small. 

Here we present first results in 4-d for the discussed method. We explic- 
itly construct an approximate solution D of (^]^) and study its properties in 
different ensembles of quenched SU(3) gauge configurations. We find that 
already with relatively few terms we obtain a good approximate solution 
for the Ginsparg- Wilson equation. We discuss properties of the spectrum 
of our Dirac operator and also analyze the role of improving the gauge field 
action. We demonstrate that improving the gauge field action considerably 
improves the spectral properties of our D. 

We would like to point out three aspects of our approximate solutions of 
the Ginsparg- Wilson equation: 



In an earlier study |15] in two dimensions we had constructed an ap- 



proximate solution of i'Ll) using the outlined ideas. It turned out that 
essential features of chiral symmetry in the lattice Schwinger model 
were properly described by our approximate solution of the Ginsparg- 
Wilson equation. It will be interesting to see if our D is able to capture 
chiral symmetry also in four dimensions. 

Our study should help to shed light on the importance of different 
terms in solutions D of the Ginsparg- Wilson equation. In particular 
for the perfect action program an economical but still rich enough 



parametrization of D is essential [16 



Another application of our D is its use as a starting operator Dq in 
Neuberger's construction. Using an approximate solution of ( |1 . 1[ ) in 



the projection (L2) has been proposed previously ||Tj and it is believed 
that an improved Dq speeds up the evaluation of the square root and 
helps to overcome the problems with the potential singularity of ( |1.2| ) 
for small eigenvalues of A. 

The article is organized as follows: In Section 2 we briefiy discuss the ex- 
pansion of D and derive the system of equations equivalent to the Ginsparg- 
Wilson equation. We discuss in Section 3 the boundary conditions and find 
a solution of the system of coupled equations. In Section 4 we analyze the 
spectral properties of our approximate solution D for ensembles of quenched 
gauge field configurations. We start with the free case in Section 4.1 followed 
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by gauge field ensembles generated with the standard Wilson gauge action 
(Section 4.2). Finally in Section 4.3 we investigate how improving the gauge 
field action can further improve the chiral properties of the lattice Dirac 
operator. The article closes with a discussion in Section 5. In particular we 
will come back to the above mentioned three possible applications of our 
Dirac operator. 



2 Constructing approximate solutions of the 
Ginsparg- Wilson equation 

To make the article self contained we briefly repeat the basic idea of our 
construction presented in In a first step we discuss the systematic 

expansion of the most general Euclidean Dirac operator in a series of simple 
operators on the lattice. Subsequently we derive the system of quadratic 
equations for the expansion coefficients which is equivalent to the Ginsparg- 
Wilson equation. 



2.1 Expansion of the most general D 

Typically the derivative term on the lattice is discretized by the following 
nearest neighbor term (we set the lattice spacing to 1): 



1 



XI 7m U^{x)6:,+f,^y - Uf,{x - fi) Ox-fL,y 



(2.1) 



However, it is perfectly compatible with all the symmetries to instead dis- 
cretize the derivative term using 



1 ^ 



7m 



Ufj,{x)U^{x+ij,) 5.j;+2fj.,y - U^{x-fi) ^Ufj.{x-2fi) ^5, 



x-2il,y 



(2.2) 

and there are many more terms one could think of. Thus an ansatz for the 
most general Dirac operator D must allow for a superposition of all of the 
possible discretizations for the derivative term. 

The terms in and ( |2.2[) can be characterized by simple paths on the 



lattice. The first example (2.1) consists of a single hop in positive /U-direction 
with a plus sign and a single hop in negative direction with a minus sign. 
Similarly the second example consists of two hops in positive (negative) //- 
direction. We now introduce a shorthand notation for such paths using an 
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ordered list of the directions of the hnks of the path. We denote a path of 
length n on the lattice by 



<ll,l2,...ln>, (2.3) 

with the li giving the directions of the subsequent hops, i.e. /j G {ilj i2, ±3, 
±4}. It is implicitly understood, that each link is dressed with the corre- 
sponding link variable U^{x). Using this notation the terms from the two 
examples (2.1) and ( |2.2| ) are denoted as 



^Etm E ^(0 <^>, (2-4) 

M Z=±At 



and 



tE^ E ^(0 <hl> ■ (2-5) 

We use the abbreviation s{l) for sign(/). Due to translation invariance the 
form of the derivative terms, and thus of the paths, is the same at all lattice 
points such that in our notation no reference to the starting point for the 
paths is necessary. 

In order to remove the doublers one also needs a term which in momen- 
tum space can distinguish between = Q (physical modes) and Pfj, = tt 
(doublers). Such a term is provided by the standard Wilson term. Due to 
the symmetries this term has to come with I in spinor space. Again we will 
allow for all possible terms. We generalize our D further, by including all 
terms also for the remaining elements Fq of the Clifford algebra, i.e. ten- 
sors, pseudovectors and the pseudoscalar. Thus the emerging lattice Dirac 
operator has the following form: 

16 

^ = Er« E <h,i2,.:iipi> ■ (2.6) 

a=l p6-P" 

To each generator r„ of the Clifford algebra we assign a set of paths p, 



each p given by some ordered set of links < Zi, •••^|p| > (compare (2.3)) 
where \p\ denotes the length of the path p. Each path is weighted with some 
complex weight Cp. 

The next step is to impose on D the symmetries which we want to 
maintain: Translation and rotation invariance and invariance under C and P. 
In addition we require our D to be 75-hermitian, i.e. we require 75-D75 = . 
This property can be seen to correspond to what leads to the CPT theorem 
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in Minkowski space, i.e. the vector generators 7^ come with a derivative 
term etc. 

Translation invariance has already been briefly mentioned above and 
requires the sets "P" of paths and their coefficients to be independent of the 
starting point. Rotation invariance implies that a path and its rotated image 
have the same weight. Parity implies that for each path p (with coefficient 
Cp) we must include the parity-reflected copy with coefficient Sp^^.^^y ■ Cp 
where the signs Spa^jty are defined by 74rQ,74 = Sp^^^fy ■ Fa- 

Of importance are the symmetries C and 75-hermiticity. It is easy to see 
that both of them imply a relation between the coefficient for a path p and 
the coefficient of the inverse path p~^. Implementing both these symmetries 
restricts all coefficients Cp to be either real or purely imaginary. Furthermore 
we find that the coefficient for a path p and the coefficient for its inverse 
p~^ are equal up to a sign Sc/iarge defined by CTaC = Scharge ' where T 
denotes transposition and C is the charge conjugation matrix. 

When implementing all these symmetries we find that paths in our ansatz 
become grouped together where - up to sign factors - all paths in a group 
come with the same coefficient. We can now write down our most general 
Dirac operator on the lattice in the form (compare ||l^ and see also [16| for 
an equivalent derivation using a slightly different notation): 

D =l[si<> +S2Y<h> <h,l2> + S4j2<h,h> ••• 

+ ^^1< h> +V2 J2 ^< ^i'^2 > + < hJi >] 

M h=±tJ. i27^±Ai 

+ V3< li,li > ... 



IJ,<1/ il=±M 



+ Y ^f^^^^P Y S{l2) S{k) Y ^ijk ai<li,lj,lk> 

IJ.<U<p ii=±M,!2 = ±'' i,j,k=l 



i,j,k,n=l 



!i=±i,;2=±2 



(2.7) 



By e we denote the totally anti-symmetric tensors with 2,3 and 4 indices. We 
choose the normalization of the elements of the Clifford algebra such that 
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the elements appear as all possible products of the 7^ without any extra 
factors of i. For this normalization the symmetries C and 75-hermiticity 
render all coefficients Si,Vi,ti,ai and pi real. To be specific, we use the 
Euclidean chiral representation for the 7^ = 7^^. In principle it would be 
possible to generalize D further by multiplying each term with a polynomial 
of traces of gauge field variables around closed loops on the lattice, but we 
do not include this possibility here. 

The above mentioned structure of paths appearing in groups is obvious. 
The paths in each group are related by symmetries and up to the sign factors 
have to come with the same real coefficient. All paths within a group have 
the same length. 

It has to be stressed, that in ( p.7| ) for each generator we show only the 
leading terms of an infinite series of groups of paths. The dots indicate that 
we omitted groups with paths that are longer than the terms we display. 
It is known |18| that no ultra-local solutions of the Ginsparg- Wilson equa- 
tion exist and thus an expansion for a solution of (|1.1D necessarily contains 
infinitely many terms. Eventually we will work with a finite D (the cor- 
responding truncation will be discussed in Section 3) but for the moment 
we keep deriving the method in its most general form, i.e. containing no 
truncation. 



2.2 The system of coupled equations corresponding to the 
Ginsparg- Wilson equation 

Let us now insert our expanded Dirac operator D into the Ginsparg- Wilson 
equation. To that purpose we multiply with 75 from the left, bring the 
terms linear in D to the right-hand side and define: 

E = - D - 75D75 + 75^75^^ • (2.8) 

We remark that E is hermitian since we implemented 75-hermiticity for 
D. Finding a solution D of the Ginsparg- Wilson equation corresponds to 
having E = 0. Evaluating the linear part of E is straightforward: When 
evaluating the product 75-D75 we find that the terms with an odd number 
of 7^, i.e. vector- and pseudovector terms pick up a minus sign, while the 
other terms remain unchanged. Thus when adding the two linear terms we 
find that the terms with an odd number of 7^ cancel, while the other terms 
pick up a factor of 2. 

The next step is to compute the quadratic term j^Dj^D. Here we have 
to multiply the various terms appearing in D. Each term is made out of 
two parts, a generator of the Clifford algebra and a group of paths. The 
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multiplication of two of these terms proceeds in two steps: First the two 
elements of the Clifford algebra are multiplied giving again an element of 
the algebra. In the second step we have to multiply the paths of our two 
terms. This multiplication can be noted very conveniently in our notation, 
where multiplication of two paths simply consists of writing the paths into 
one long path: 

<h,l2---ln> X <l'l,l2---ln' > = < h^2--ln,l'lA--A' > ■ (2-9) 

It is straightforward to establish this rule by translating back to the algebraic 



expression of our examples ( |2.lD ,(2.2) and performing the multiplication in 
this notation. It can happen that after multiplying two paths a hop in some 
direction li is immediately followed by its inverse —li. These two hops then 
cancel each other and we find 



< h...l 



1 • ••'i— 1 ) 



^7 1 1 ^7 



i 1 • 



■In > 



< U...I 



■In > 



(2.10) 



This rule is used to reduce all products of paths appearing in 'j^D^^D to 
their true length. In a final step we decompose the product terms into 
groups related by the symmetries in the same way as we did above when 
constructing the most general ansatz for D. Adding the linear and quadratic 
terms of ( p.8[ ) we end up with the following expansion for E: 

E =l[e^<> +e^^</i> +e^^ <hM> + elY,<hA> ■■■ 

h h^h h 



h=±^l 



i,j,k=l l4yt±fl,U,p 



i3 = ±P 



— < Z4, li, Ij, Ik >] + e2[< li, li, lj,lk > ~ < Ijj Ik 

4 

+ 75 E s{li) S{l2) s{h) S{U) ^ eijkn el <li,lj,lk,ln> 

'i=±i''2=±2 i,j,k,n=l 
!3 = ±3,i4=±4 



(2.11) 



We remark that all the algebraic steps leading to the expansion of E are 
straightforward to formalize and for the higher orders in the expansion we 
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used a computer program. As for D, the expansion of E is an infinite series 
and we display here only the leading groups of paths. The coefficients ef 
are now quadratic polynomials in the original coefficients Si,Vi,ti, ai and pi 
given by 



: -2si + sl + 8sl + 48sl + 8sl + 8vj + 96v| + 8v| + 48tj + 192a^ 
+ 384p? ... , 

62 = - 2s2 + 2siS2 + I2S2S3 + 2s2S4 + 12wiW2 + 2t;it;3 ... , 

63 = - 2s3 + 2siS3 + S2 + 4s| + 2S3S4 + 4^2 + 2t;2l'3 ... , 

el = - 2s4 + 2siS4 + S2 + 6si - - 6t? -24a? - 48p? ... , 



el = 


- S2V1 - 4S3W2 


- 2S4V2 - S3V3, - vsh - Av2tl ... , 




- 2ti + 2siti - 

+ letipi ... , 


2s4ti -vf- Avl - 2V2V3 - + 8viai 


e\ = 


— S2ai + V2tl — 


V3P1 ■■■ , 


e2 — 


- V2tl ... , 




a 

es = 


- S2ai - 2v2Pi 




64 — 


- 2s2ai + 2^2*1 


- 4:V2Pi ... , 


ei — 


- 2pi + 2sipi - 


2s4pi — 2viai + ti ... . 



8a? 



(2.12) 

Due to different symmetry properties {E is hermitian, D is 75-hermitian) 
there are terms in E which do not occur in D and vice versa. Each coefficient 
itself is an infinite series of terms. For a solution of the Ginsparg- Wilson 
equation we must have E = 0. It is easy to see, that the groups of paths 
appearing in the expansion ( 2.11| ) are linearly independent and hence for 



= all coefficients ef have to vanish simultaneously. Thus we have 
rewritten the problem of finding a solution of the Ginsparg- Wilson equation 
to solving the system (2.12) of coupled quadratic equations (set all right- 
hand sides ef = 0). 



3 Solving the system of coupled equations 

In the last section we have shown that the Ginsparg- Wilson equation is 
equivalent to a system of coupled quadratic equations for the expansion 
coefficients of D. In this section we now truncate our expansion of D and 
find solutions for the coupled equations. Before we do this let us discuss the 
boundary conditions which have to be added to the system ( 2.12| ). 
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3.1 Boundary conditions 

For the free case the situation is simple: In this case we can compute the 
Fourier transform D{p) of D. For smah momenta the massless Dirac oper- 
ator should obey 

D{p) = + 0{p^) . (3.1) 

This leads to two more supplementary equations for the coefficients, one 
for the constant term to vanish and the second one sets the slope of the 
dispersion relation equal to one: 

= si + 8 S2 + 48 S3 + 8 S4 ... , (3.2) 

1 = 2ui + 24t>2 + 4v3 ... . (3.3) 



We remark, that (3.2) is implicitly guaranteed by the Ginsparg- Wilson equa- 
tion, since the Ginsparg- Wilson circle^ runs through the origin. When we 
construct approximate solutions below it will be necessary to explicitly im- 
plement (|3.2[) . When we started to solve the system of equations ( |2.12|) 
respecting the boundary conditions ( |3.2| ), ( ^^ ) we found, that for this set- 
ting one has to include a relatively large number of terms in the expansion 
( |2.7D of D. Since every new term in D drives up the cost of a numerical treat- 
ment of our Dirac operator we decided to allow for an additional freedom 
in our equations. The idea is to allow some /^-dependence by introducing a 
new parameter, similar in spirit to the variable s in the overlap construction 



(1.2), which may be used to improve locality of D , i.e. to make higher terms 



in the expansion (2.7) of D less important. At the same time we do not want 



to change our coupled equations ( |2.12|) since they guarantee that - up to 
a certain order - we deal with a solution of the Ginsparg- Wilson equation. 
Thus we only modify the boundary conditions (|3.2| ), (^). 

Before we discuss the modifications let us first develop the idea for the 
case of the standard Wilson action. Wilson's lattice Dirac operator Dyy is 
a sum of a constant and a hopping term H 

Dw = 4: - K H , 

(3.4) 



^It is well known, that an exact, 75-hermitian solution of the Ginsparg- Wilson equation 
has its spectrum on the so-called Ginsparg- Wilson circle, i.e. the circle of radius 1 with 
center 1 in the complex plane. 
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For the free case the Fourier transform of the hopping matrix is given by 
H{p) = 8 — i 2 ^ + 0{p'^). Thus the boundary condition corresponding to 



Q reads 

and from this equation one finds that for the free case the correct value of the 
parameter is given by k = 1/2. In fact, in the notation of our parametriza- 
tion of D we have si = 4, S2 = — k and vi = k with ah other coefficients 
vanishing. The condition (|3.5|) for the Wilson operator therefore corresponds 



to both our boundary conditions (|3.2D and (|3.3|) . In principle one could gen- 



eralize the Wilson action and lift the degeneracy between —S2 and vi which 



would amount to two boundary conditions similar to ( ^.21) , (3.3) instead of 
the single condition (3.5). 



The situation changes when we include gauge fields. One finds that 
in order to drive the system critical one needs to change the value of k. 
E.g. when analyzing quenched SU(3) gauge theory at /? = 6.0 one finds 
K = 0.624 (see below). Obviously this value of kappa is not a solution of the 
boundary condition (3.5) but instead solves the more general equation 



4 - 8kz + 2zKzj/ + 0(/) = i^+0{p^) . (3.6) 

Here z is a real function of /3. For quenched gauge fields at /3 = 6, we find 
e.g. z = 1.603 = 1/0.624. 



We now generalize the boundary conditions (3.2), ( |3.3| ) analogous to 
( |3.6D being a generalization of ( |3.5D . In the boundary condition we multiply 
each coefficient by a power z" of some real, /3-dependent number z. The 
exponent n is given by the number of hops of the corresponding term in our 
expansion (|2.7|). Thus e.g. the coefficient si, which is the constant term, 
does not get changed at all, S2 which is the coefficient for the single hop 
obtains a factor of z, S3 which corresponds to two hops becomes multiplied 
by z"^ et cetera. We allow for two different coefficients Zg and z„ in the scalar 
and vector sectors. Our generalized boundary conditions then read 

= si + 8 S2 -Zs + 48 S3 + 8 S4 zl ... , 

1 = 2 vi z^ + 2Av2 zl + A zl ... . (3.7) 

The two parameters Zs(/3) and z.u{(3) can be used to optimize the properties 
of D while at the same time maintaining unchanged the equations ( p.l2| ) 
which are equivalent to the Ginsparg- Wilson equation. 

We emphasize, that these two supplementary conditions are just a prac- 
tical means to effectively reduce the number of necessary terms in the Dirac 



10 



operator. In some sense they serve as a guiding principle to select more local 
actions in the huge space of possible actions satisfying the Ginsparg- Wilson 
condition. 

We can now restrict ourselves to a much smaller set of terms in the 
expansion ( p.TD of D. Once this set is chosen, Zg and will be determined 
by optimizing the properties of the spectrum of D near the origin. The 
parameter Zg is fixed by the requirement that the small eigenvalues of D fall 
on the Ginsparg- Wilson circle (i.e. m = 0); z^ can e.g. be fixed by requiring 
the slope of the vr-dispersion relation to be equal to 1. In this preliminary 
study we will simply determine z^ by optimizing the alignment of larger 
eigenvalues along the circle. 

It is important to remark, that already including only a few terms in the 
expansion ( ^.7| ) of D and solving the corresponding system ( |2.12 ) consid- 



erably orders the small eigenvalues of D (see below). Thus no complicated 
fine-tuning procedure is necessary for Zs and z^ and they can simply be de- 
termined by analyzing the spectrum of D on a few background gauge field 
configurations. 

Finally we remark, that another feature which one would like to imple- 
ment is 0{a) improvement It is known, that an exact solution of the 
Ginsparg- Wilson equation is already 0{a) improved. However, here we will 
discuss only an approximate solution and it is useful to require 0{a) improve- 
ment independently of the Ginsparg- Wilson equation. 0{a) improvement 
may be achieved by adding to D the so-called clover- leaf term Cgyj-^ a-F pO| ]. 



At tree level the coefficient Cgw equals the factor in front of the Laplace-type 
contribution to D. For non-perturbative improvement Csw can e.g. be de- 
termined using the Schrodinger functional [^]. Since we do not attempt to 
determine Csw non-perturbatively here we quote the condition at tree level 
(i.e. Csw = !)• It is obtained by expanding the terms in the tensor sector 
of our D for small lattice spacing, extracting their ^ cr • content and set- 
ting the factor in front of this term equal to the factor of the Laplace-term 
contribution of D. One finds (compare also 

S2 + 12 S3 + 4 S4 ... = 4 ti + 32 t2 + 16 U ... . (3.8) 



It is straightforward to include Cgw as a free parameter in (|3.8D but as already 
remarked above we do not attempt a nonperturbative evaluation of Csw here 
and thus work with the tree level equation ( |3.8D throughout this paper. 

At this point we would like to comment on the fate of the doubler modes. 
When analyzing our solutions for the free case we always find that the eigen- 
values for the doublers (at least one component of the four-momentum 
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equals vr) are located near 2 in the complex plane as is the case for exact, 
doubler-free solutions of the Ginsparg- Wilson equation. Thus our approx- 
imation scheme already takes care of the doublers and we do not need to 
supplement the boundary condition (|3.lD by an additional equation for the 
doublers. 



3.2 Truncation of D and numerical solution of the system of 
coupled equations 



In the last section the set of equations ( 2.12| ) together with the boundary 



conditions ( ^.Tj) and (|3.8|) was discussed. In principle the system ( ^.12 ) 



contains infinitely many equations, each of them with infinitely many terms. 
Also the boundary conditions contain infinitely many terms. The next step 



thus will be a truncation of the expansion (2.7) for D which will reduce the 
infinite problem to a finite problem. 

Any reasonable truncation should have a parameter which controls the 
size of the remainder of the approximation. In our case such a parameter is 
given by the length of the paths in each term in the expansion ( |2.7| ) of D. 
We find that roughly the size of the coefficients sf,vf,tf,af,pf decreases 
exponentially as the length of the paths in the corresponding terms increases 
(compare Table ^ in the appendix) . Thus neglecting terms with longer paths 
provides a natural cutoff scheme for the expansion of D. The exponential 
decrease of the coefficients not only justifies our cutoff, but is also important 
for the physics described by D. A solution of the Ginsparg- Wilson equation 
cannot be ultra- local jl^. However, in order to remain in the correct uni- 
versality class, D has to be local, i.e. Dx^y has to decrease exponentially as 
the distance |x — ?/| on the lattice increases. In Fig. || we plot the size of the 
coefficients as a function of |x — y| on a logarithmic scale. We show the coef- 
ficients for the free case (full circles are used to represent all the coefficients) 
as well as those used for the quenched gauge field ensemble (triangles) from 
the Wilson gauge action at /? = 6.0 (cf. Table |3|). The coefficients for the 
other three ensembles show a similar behavior. The size of the coefficients 
is consistent with being bound by an exponential decay with |a; — ?/|. Thus 
we have verified, that we approximate a healthy, local solution of ( |1.1D . 

Let us now come to the more practical aspects of the truncation. It 
is clear that adding more terms in the expansion ( |2.7|) of D improves the 
quality of the approximation of a solution of the Ginsparg- Wilson equation. 
On the other hand each new term drives up the numerical cost. Here we 
present an operator which manages to give a decent approximation of a 
solution of (|l.lD but at the same time is still relatively cheap to simulate. 
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Figure 1: Size of the coefficients as a function of the distance |x — y| on 
the lattice. We use fuh circles for ah the coefficients in the free case, the 
triangles represent the coefficients of the quenched ensemble aX (3 = 6.0 (cf. 
Table 0). 



When constructing the operator we started with a simple parametrization 
and added only terms which gave rise to a considerable improvement of the 
spectral properties. We stopped adding terms when our D had a satisfactory 
balance of good chiral properties and low numerical cost. 

Also the choice of equations from the set ( 2.121 ) allows for some freedom 
(the boundary conditions ( |3.7D and ( |3.8D are always implemented). This 
freedom again goes back to the non-existence [^] of ultra-local solutions 
If we could solve all equations (2.12) for a finite parametrization 
of D this would, however, amount to an ultra-local solution. The 
loophole out of this dilemma is the fact that the system ( p. 12 ) is always 
overdetermined. Thus we can only solve the equations corresponding to the 
leading terms of the expansion (2.11) of where again the length of the 
paths is the expansion parameter. 

The solution which we present here has altogether 17 terms in the ex- 
pansion (^]^) of D. The maximum path length is 4 and we only included 
terms in the scalar, vector and tensor sectors. Except for one term (sg) all 
terms have paths on the hypercube. For the vector and tensor sectors we 
allowed only terms up to length 3 on the hypercube, since longer terms in 
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these sectors create many new entries in the fermion matrix which quickly 
increase the numerical cost. A detailed description of the terms used in our 
D and the values of the coefficients are given in the appendix in Table |2| and 
Table |3|. Besides the three boundary conditions (^), (|3.8|) we implement 
14 equations from the system ( p.l2 ) to match the number of expansion co- 
efficients (17). These equations correspond to the shortest terms with paths 
on the hypercube contributing to the expansion ( p.ll[ ) of E. We remark, 
that it was straightforward to exactly solve the equations ( 2.12[ ), (^), ( |3.8| ) 
with a standard solver |22] and we also found that the system is very stable, 
i.e. we found only a single solution for our setting. 



4 Properties of our Dirac- operator 

In this section we analyze the properties of our Dirac operator for different 
ensembles of background gauge fields. We start with the analysis of the free 
case, followed by a study using quenched SU(3) configurations generated 
with the standard Wilson action for the gauge fields. We conclude with 
analyzing gauge fields from the improved Liischer-Weisz action. 



4.1 The free case 

For the free case it is possible to completely diagonalize the Dirac operator 
using Fourier transformation. In Fig. ^ we show the free spectrum of our D 
in the complex plane. In addition to the eigenvalues (symbols) we also show 
the Ginsparg- Wilson circle (full curve) which supports the spectrum of any 



exact solution of (1.1). Throughout this article we present our spectra for 
8^ lattices. This is a relatively small volume for usual standards, however, 
the Dirac operator studied here is ultralocal and its extent is well within the 
lattice size. 

Let us briefly discuss some features of the free spectrum. It is obvious, 
that our approach optimizes the alignment of the eigenvalues along the circle 
for the eigenvalues in the physical branch, i.e. the eigenvalues in the vicinity 
of the origin. The alignment of the eigenvalues in the doubler branches 
is less perfect. The important feature, however, is the clear separation of 
the doubler modes and the physical eigenvalues. We furthermore found, 
that adding additional terms in (^]^) systematically improves the situation 
also for the doublers. This property of our approach of first optimizing 
the physical branch and aligning the doubler modes along the circle when 
adding higher terms persists also when analyzing the spectrum for non- 



trivial background gauge fields and was observed also for the 2-d case |15]. 
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0.0 0.5 1.0 1.5 2.0 



Figure 2: Spectrum of our Dirac operator for the free case in the complex 
plane. The represent the eigenvalues. We also show the Ginsparg- Wilson 
circle, i.e. the circle of radius 1 and center 1. 



At this point it is interesting to compare our method to the study of 
truncated perfect actions for free fermions |17, For free fermions the 

perfect action can be computed explicitly Pfas an infinite series and is a so- 
lution of the Ginsparg- Wilson equation. After truncation to e.g. only terms 
on the hypercube the result is an approximate solution of the Ginsparg- 
Wilson equation for the free case. 

This truncated perfect action will however differ from our approach in 
two important aspects: Firstly actions for free fermions only contain scalar 
and vector terms. Terms in the tensor, pseudovector and pseudoscalar sec- 
tors vanish identically when no gauge field is coupled. Thus e.g. terms in 
the tensor sector which are necessary for 0{a) improvement and also play 
an important role for a smooth spectrum in the physical branch |25| have to 
be included a posteriori. Secondly for each possible endpoint of a fermion 
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path the truncated perfect action for free fermions provides only a single co- 



efficient. On the other hand, our symmetry analysis (compare (2^)) shows, 
that different paths leading to this endpoint can come with different coeffi- 
cients. Thus the distribution of the single coefficient of the free perfect action 
among the various coefficients allowed by the exact symmetry analysis is an 
additional task. 

Overcoming these problems, Bietenholz [|2^ presented interesting results 
for an approximate solution of the Ginsparg- Wilson equation in 4-d. After 
introducing fat links and a link amplification factor for the coefficients of the 
paths, an approximate solution of the Ginsparg- Wilson equation also for the 
case with gauge fields is obtained from the truncated perfect action. It is 
interesting to note that the spectra presented in |24] show a behavior quite 



contrary to the spectra for our D: They come with a very smooth doubler 
branch while the modes in the physical branch have quite large fiuctuations. 
Bietenholz expresses confidence that the physical branch can be smoothened 
by adding tensor terms but expects that these terms will at the same time 



destroy the good alignment of the eigenvalues in the doubler branch | . 



4.2 Results for quenched gauge configurations generated with 
the Wilson gauge action 

After having studied the properties of our D for the free case we are now 
analyzing the Dirac operator in ensembles of quenched SU(3) gauge field 
configurations. We concentrate on the branch of the spectrum near small 
eigenvalues, which is the one most relevant in the continuum limit. We 
compute the eigenvalues in this physical branch of the spectrum and compare 



them with the eigenvalues for the Wilson operator (3.4). The values of k 
used for the Wilson operator (^^) can be found in Table |^ and will be 
discussed below. 

For the case of non-trivial background gauge fields the Dirac operator 
can no longer be diagonalized completely unless one uses very small lat- 



tices. Here we use the Implicitly Restarted Arnoldi Method [26| to compute 
eigenvalues in the physical branch of the spectrum, i.e. near the origin. The 
method allows to specify a search criterion for the eigenvalues and we use this 
feature to compute the eigenvalues with the smallest real parts. In order to 
get insight into the typical behavior of the spectrum of our D, we analyzed 
the eigenvalues on 4 ensembles of quenched gauge fields. Each ensemble 
consists of 20 well decorrelated configurations. Two of the ensembles were 
generated with the standard Wilson action at /3 = 6.0 and (3 = 5.85. They 
will be discussed in this subsection. The other two ensembles were generated 
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Figure 3: The physical branch of spectra of Dirac operators in the complex 
plane for quenched SU(3) gauge field configurations generated using the 
Wilson gauge action. From left to right we show: The spectrum of our D 
at P = 6.0, for Wilson's Dirac operator at /? = 6.0, our D at (3 = 5.85 and 
finally for Wilson's D at (3 = 5.85. 



using the improved Liischer-Weisz action |27] and will be analyzed in the 
next subsection. For each gauge field configuration we computed 60 eigen- 
values in the physical branch for both our Dirac operator as well as for the 
Wilson-Dirac operator. For the plots in Fig. ^ and Fig. ^ we increased the 
number of computed eigenvalues to 200. All these calculations were done on 
8^ lattices with anti-periodic boundary conditions for the Dirac operator. 

In Fig. ^ we show the physical branch of the spectrum of our operator 
for quenched gauge field configurations generated with the standard Wilson 
gauge action at /3 = 6.0 and (3 = 5.85. We compare these spectra with 
the spectra of the Wilson-Dirac operator on the same configurations. The 
symbols are the numerically computed eigenvalues and the full curve is the 
Ginsparg- Wilson circle. 

It is obvious, that for our D the eigenvalues are more ordered than for 
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Wilson's D for both /3 = 6.0 and /? = 5.85. The spread of the eigenvalues 
of our D is reduced when compared to the Wilson case and the eigenvalues 
tend to order near a curve. However, what is also obvious is the fact that 
although this curve touches the origin, it also bends away from the Ginsparg- 
Wilson circle as the size of the imaginary parts increases. As can be seen 
from the Fourier transform in the free case, it is essentially the vector terms 



in the expansion (2^) contributing to the imaginary parts of the eigenvalues. 
As we have already remarked above, we allowed only for a relatively small 
number of vector terms (compare Table in order to reduce the numerical 
cost. We thus attribute the bending away of the eigenvalues to this reduced 
approximation of the vector sector. In the next section we will, however, 
demonstrate, that our parametrization of the vector sector is sufficient if one 
uses improved gauge field actions. 

It is also interesting to study the behavior of the real eigenmodes - the 
P = 5.85 part of Fig. ^ shows one real eigenvalue. Typically the fluctuations 
of the real eigenvalues of the Wilson-Dirac operator are relatively large. In 
particular for so-called exceptional configurations (in the sense of a break- 
down of the matrix inversion) a real mode has fluctuated so heavily towards 
small values that it compensates for the bare quark mass and the resulting 
zero mode leads to a breakdown of the matrix inversion. We found, and 
Fig. ^ shows an example, that our D strongly suppresses the fluctuations 
of the real modes and turns them into almost perfect zero modes. We thus 
expect that at nonzero bare quark mass our D considerably reduces the 
numerical difficulties with exceptional configurations. 

Let us now discuss the overall fluctuation of the rest of the physical 
branch of the spectrum. An exact solution of the Ginsparg- Wilson equation 
has its spectrum on the Ginsparg- Wilson circle and only the density of 
eigenvalues on this circle will fluctuate. The eigenvalues of the Wilson- 
Dirac operator, on the other hand, have relatively large fluctuations in all 
directions. When analyzing the ensemble of 20 configurations for both P = 
6.0 and /? = 5.85 we found that for our D the fluctuations of the physical 
edge of the spectrum are considerably suppressed when compared to the 
Wilson-Dirac operator. In order to quantify this statement we fit a circle of 
the form 



y = ±i^Jl - (1 + e - xf , (4.1) 

to the 10 smallest eigenvalues in the physical branch omitting the exactly 
real eigenvalues, because they deviate significantly in the case of the Wilson 



operator (compare above). Eq. ( [O] ) describes a circle of radius 1 around 



(1 + e, 0) in the complex plane. The parameter e of the fit gives the position 
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where the circle crosses the real axis. The fluctuation of e provides a measure 
for the fluctuation of the eigenvalues in the physical branch of the spectrum. 



|3 = 6.0 P = 5.85 




8 8 

Figure 4: Histograms indicating the fluctuation of e (i.e. the fluctuation of 
the physical branch of the spectrum) for two ensembles of 20 quenched gauge 
field configurations generated with the Wilson gauge action at /? = 6.0 (left 
column) and (3 = 5.85 (right column). The top row shows the results for 
our D, the bottom row displays the results for the Wilson-Dirac operator. 

In Fig. ^ we show for each of our ensembles a histogram of the values 
of £ for all 20 configurations, comparing the results for our Dirac operator 
with those for Wilson's D. It can be seen, that our D has considerably 
reduced the fluctuations of the physical edge of the spectrum, in particular 
for /? = 6.0. 

At this point we can also comment on the determination of the critical 
K for the Wilson-Dirac operator (|3.4| ). Since the first term is proportional 
to the unit matrix, it is sufficient to diagonalize the hopping matrix H. The 
value of K, was then adjusted such that the distribution of e is centered at 
0. The values of k obtained for our 4 ensembles are given in Table |l] below. 
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free case Wilson action Liischer-Weisz action 

(3 = 6.0 13 = 5.85 (5i = 8.45 j3i = 8.15 

i < C/pZ > O 0.594(1) 0.576(1) 0.652(1) 0.633(1) 

K 0.5 0.624 0.642 0.607 0.623 



Table 1: Parameters for our ensembles of quenched gauge field configurations 
and the free case. We list the expectation value < Upi > /3 of the plaquette 
and the critical k we use for the Wilson operator. 

4.3 The effect of improving the gauge action 

In this section we discuss the effect of improving the gauge field on the 
spectral properties of the Dirac operator. In the context of domain wall 
fermions it has been found that improving the gauge field action allows 
to work with a smaller extension of the 5-th direction ||2^. Also for the 
4-d setting it is worth testing if improving the gauge field action leads to 
an improvement of the spectral properties of D in particular since from 
a numerical point of view improving the gauge action is an inexpensive 
measure. 

In order to analyze the effects of gauge improvement we generated 2 
ensembles of 20 gauge configurations each, using the Liischer-Weisz action 
[p7| , [29| ]. The parameters were adjusted such that the physical scale (the 
effective string tension) approximates the two scales of our ensembles from 
the standard Wilson gauge action at /? = 6.0 and (3 = 5.85. Our estimator 
for the effective string tension is based on Wilson loops up to extent 3 only 
and therefore not to be compared with the values derived on larger lattices 
with a detailed finite size analysis. However, our principal incentive was 
not to obtain such a precise value, but to identify the values of the gauge 
couplings where the two gauge actions have roughly the same scale. 

To be specific, we use the setting for the improved gauge action as pre- 
sented in [29 1 but with only the rectangle term. Explicitly, the gauge field 



action reads 

S[U] = PiYl Tr(l - Upi) + /?2 E l^"" Tr(l - Urt) , (4.2) 

pi rt 

where the first sum is over all plaquettes and the second sum over all 2 x 1 
rectangles. (3i is the driving parameter while P2 can be computed from /?i 
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using tadpole improved perturbation theory |3C] giving [2£] 



(^2 = - [1 + 0.4805 a,] . (4.3) 

20 u 







with 



1^ ^ v\V4 



,x X.,, In (iRe Tr C/p, ) 

wo = i-ReTr{U„i)) ' , a, = (4.4) 

V3 \ P'V ' 3.06839 ^ ^ 

The coupling (3i is determined self-consistently with uq and for a given /3i . 
We adjust the parameter Pi such that the physical scale approximates the 
two scales of our ensembles from the standard Wilson gauge action at /? = 6.0 
and P = 5.85, i.e. such that the effective string tension of the improved 
ensemble (as discussed above) roughly matches the value of the ensemble 
from the standard Wilson action. The resulting values are Pi = 8.45 which 
corresponds to P = 6.0 for the Wilson action and Pi = 8.15 corresponding 
to P = 5.85. We list our results for the plaquette expectation value and the 
values of the critical k determined for the Wilson-Dirac operator in Table ||. 

Like for the ensembles generated with the Wilson gauge action also here 
we computed for each configuration 60 eigenvalues in the physical branch for 
both, our D as well as for Wilson's Dirac operator. When inspecting these 
eigenvalues, one finds that the spectral properties of our D are considerably 
improved. In Fig. ^ we show examples of spectra for one configuration from 
each of the ensembles. As in the last section we show 200 eigenvalues in the 
physical branch for our D as well as for the Wilson-Dirac operator both for 
the same gauge configuration. 

When comparing these plots to Fig. ^ one finds that the eigenvalues of 
our D are somewhat better aligned along a single curve, i.e. the distribution 
is narrower. Also the imaginary parts are larger and the curve is following 
the Ginsparg- Wilson circle much closer than it was the case for the stan- 
dard Wilson gauge action. The improved spectral properties of D can be 
attributed to the following feature of the improved gauge field ensemble: 
When comparing the plaquette expectation value of the improved ensemble 
with that of the Wilson ensemble at the same physical scale one finds that 
improvement brings the plaquettes considerably closer to 1 (see Table |^) . A 
second, although smaller effect, might be the better rotational properties of 
the gauge fields from the improved action p^ . 

We remark, that also the eigenvalues of the Wilson-Dirac operator are 
slightly more ordered when using the improved gauge field action. It is 
expected that this already could help to improve on the problems with the 



singularities in the inverse square root of the overlap operator |31|. 
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Figure 5: The physical branch of spectra of Dirac operators in the complex 
plane for quenched SU(3) gauge field configurations generated using the 
Liischer-Weisz gauge action. From left to right we show: The spectrum of 
our D at (3i = 8.45, for Wilson's Dirac operator at /?i = 8.45, our D at 
/?! = 8.15 and finally again Wilson's D at Pi = 8.15. 



We conclude this subsection with repeating last section's analysis of the 
fluctuations of the physical branch, now for the improved ensemble. When 
comparing Fig. ^ and Fig. Q one finds, that for both our D as well as for the 
Wilson-Dirac operator the distribution of e becomes slightly narrower when 
using the improved gauge action. 

We thus find that using the improved gauge action is a numerically in- 
expensive measure leading to a considerable improvement of the spectral 
properties of the lattice Dirac operator. In particular we observe that im- 
proving the gauge fields allows to work with a relatively small number of 
terms in the vector sector of D and still obtain a good approximation of a 
solution of the Ginsparg- Wilson equation. 



22 



|3^ = 8.45 P^ = 8.15 




Figure 6: Fluctuation of the physical branch of the spectrum for the two 
ensembles from the Liischer-Weisz action at Pi = 8.45 (left column) and 
/?i = 8.15 (right column). The top row shows the results for our D, the 
bottom row displays the results for the Wilson-Dirac operator. 



5 Discussion 

In this article we have presented first tests in 4-d for a new method of 
constructing approximate solutions of the Ginsparg- Wilson equation. The 
most general Dirac operator on the lattice is systematically expanded in a 
series of simple operators. The Ginsparg- Wilson equation turns into a set 
of coupled quadratic equations for the expansion coefficients. For a finite 
parametrization of D these equations can be solved and the corresponding 
D is an approximate solution of the Ginsparg- Wilson equation. We im- 
plement boundary conditions which allow to work with a very economical 
parametrization of D and still obtain a good approximation of a solution 
of the Ginsparg- Wilson equation. Our framework allows to systematically 
include 0{a) improvement. 

We find that our method has the intriguing feature that already a few 
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terms lead to a good alignment on the Ginsparg- Wilson circle for the phys- 
ical modes. The alignment of the doubler branch can be systematically 
improved by adding additional terms in the parametrization of D. We have 
demonstrated that the fluctuations of the eigenvalues in the physical branch 
are much smaller than e.g. for the Wilson-Dirac operator and we expect that 
our D gets essentially rid of the problems with exceptional configurations. 
Finally we have observed that using an improved action for the gauge fields 
considerably improves the spectral properties of our D. 

At this point we would like to comment on the numerical cost of our 
Dirac operator. So far we have not implemented an optimized matrix-vector 
multiplication for our D. Our test program was kept general in order to be 
able to test Dirac operators with more terms than the ones we present here. 
Thus we give a theoretical estimate of the numerical cost of our Dirac op- 
erator and compare it to the cost of the standard Wilson-Dirac operator. 
To leading order the numerical cost is determined by the number of non- 
vanishing entries in the fermion matrix. For our D this number is given by 
297 where we have factored out the volume and a factor of 12 for color and 
spinor structure. The corresponding number for Wilson's Dirac operator is 
17. Thus to leading order we find a factor of 297/17 = 17.47 when comparing 
the numerical cost for our D with Wilson's Dirac operator. This estimate, 
however, assumes that the all necessary products of the link variables can 
be stored. This is possible for smaller lattices up to 10^, but for larger lat- 
tices only pieces of paths can be stored. With a good storage strategy the 
additional cost for building up the longer paths can be kept below a factor 
of 2. We estimate our D to be approximately 25 times more expensive when 
compared to Wilson's D. Finally, wc also found that the Arnoldi diagonal- 
ization routine we use typically needs 35 % less matrix multiplications for 
our D as compared to Wilson's Dirac operator. 

Let us now come back to the three possible applications of our D which 
we mentioned in the introduction: 

• We think that it would be interesting to use our D for analyzing the 

pion spectrum in the quenched approximation. The fact that the fluc- 
tuations of the spectrum near the origin are highly suppressed should 
allow one to perform simulations at small bare quark masses without 
running into problems with exceptional configurations. Also studies of 
the distribution of the eigenvalues of our D could be used to compute 
the chiral condensate from the Banks-Casher formula. 

• Our successful construction of an approximate solution for the Ginsparg- 
Wilson equation also seems to be good news for the project of con- 
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structing perfect fermion actions using block spin transformations. It 
seems quite feasible to construct a sufficiently rich parametrization 
with only 15 to 20 parameters. It seems to us that in the scalar sector 
it is necessary to include relatively many terms, while in the higher 
sectors one can try to be more economical, in particular when using 
improved or perfect gauge actions. 

• Finally our D is also a good candidate for an improved starting op- 
erator Dq in the overlap projection (p..2|): It is already much closer 
to a solution of the Ginsparg- Wilson equation than e.g. the standard 
Wilson-Dirac operator which is usually used for Dq. Thus we expect 
that an expansion for the inverse square root converges faster when 
using our D as starting operator. 
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A Technical appendix 



In this appendix we describe in more detail the terms in our Dirac operator 
and give the values for the coefficients which we use for the 4 ensembles of 
quenched gauge field configurations. 

As has been pointed out in Section 2.1, the most general Dirac operator 
D can be expanded in the series (|2.7| ). Each term in this series is character- 
ized by three pieces: 

1. a generator of the Clifford algebra, 

2. a group of paths, 

3. a real coefficient. 

The paths within a group can have different signs which are determined by 
the symmetries, C, P, 75-hermiticity and rotation invariance (for the oper- 
ation of these symmetries see Section 2.1). The symmetries also determine 
which paths are grouped together. Thus it is sufficient to characterize a 
group of paths by a single generating path and all the other paths in the 
group as well as their relative their sign factors can be determined by ap- 
plying the symmetries. 

In addition, for the vector and tensor terms appearing in our D it is 
sufficient to give the paths only for one vector (tensor) since rotation invari- 
ance immediately fixes the structure for the other vector (tensor) terms. In 
order to describe our D we start with listing the three determining pieces 
for each term in Table ^. 

In Section 4 we analyze our D in different ensembles of quenched back- 
ground gauge fields. In particular we use ensembles generated with the 
standard Wilson action at /3 = 6.0 and /? = 5.85 and ensembles from the 
Liischer-Weisz action at (3i = 8.45 and /3i = 8.15 (for more details see Sec- 
tion 4.3). For all these configurations we use the same parametrization of D. 
Only the values of the coefficients sf,vf and tf and the two normalization 
factors Zg and Zy for the boundary conditions differ. In Table |3| we list their 
values. 
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Clifford generator Generating path Name of coefficient 



I 


<> 


Sl 


I 


< 1 > 


S2 


I 


< 1,2 > 


S3 


I 


< 1,2,3 > 


S5 


I 


< 1,1,2 > 


se 


I 


< 1,2,-1 > 




I 


< 1,2,3,4 > 


sio 


I 


< 1,2,-1,3 > 


Sll 


I 


< 1,2,-1,-2 > 


Sl3 


71 


< 1 > 


Vl 


71 


< 1,2 > 


V2 


71 


< 1,2,3 > 


V4 


71 


< 2,1,3 > 


V5 


71 72 


< 1,2 > 


h 


71 72 


< 1,2,3 > 


t2 


71 72 


< 1,3,2 > 


t3 


7i72 


< 1,2,-1 > 


t5 



Table 2: Description of the terms in our D. 
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free case Wilson action Liischer-Weisz action 







(3 = 6.0 


13 = 5.85 


f3i = 8.45 


f3i = 8.15 


Sl 


+1.488513 


+ 1.541745 


+ 1.536335 


+1.545142 


+1.546343 


S2 


-0.030753 


-0.066240 


-0.069085 


-0.061723 


-0.063831 


S3 


-0.011132 


-0.013808 


-0.014036 


-0.014045 


-0.013612 


S5 


-0.002128 


-0.002763 


-0.002844 


-0.002571 


-0.002690 


se 


-0.000691 


+0.002396 


+0.002649 


+0.002228 


+0.002198 


S8 


-0.005842 


-0.005813 


-0.005786 


-0.005415 


-0.005811 


SlO 


-0.000303 


-0.000557 


-0.000593 


-0.000512 


-0.000526 


Sll 


-0.000830 


-0.001217 


-0.001241 


-0.001181 


-0.001200 


Sl3 


-0.010061 


+0.007853 


+0.007918 


+0.007831 


+0.007788 


Vl 


+0.032416 


+0.091060 


+0.096935 


+0.107048 


+0.087543 


V2 


+0.032158 


+0.022333 


+0.022054 


+0.017838 


+0.022207 


VA 


+0.008365 


+0.006828 


+0.006007 


+0.007926 


+0.007540 


Vb 


-0.013326 


+0.000299 


+0.000423 


+0.001779 


-0.000364 


h 


-0.108727 


-0.103545 


-0.103318 


-0.102314 


-0.103576 




-0.002129 


-0.003544 


-0.003715 


-0.003206 


-0.003394 


h 


+0.001356 


+0.002333 


+0.002444 


-0.002198 


+0.002239 




-0.006355 


-0.006945 


-0.007266 


-0.006170 


-0.006645 




1.0 


0.876 


0.866 


0.895 


0.885 


Zy 


1.0 


0.868 


0.885 


0.846 


0.854 



Table 3: The numerical values of the coefficients sf , vf and tf and the factors 
Zg and Zv for the different ensembles of gauge fields. 
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